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The  principal  investigators,  Greg  Forest  and  Stephen  Bechtel,  have  made  progress 
on  analytical,  numerical,  and  experimental  fronts  in  our  studies  of  non-Newtonian  fluids. 
The  resources  afforded  by  this  grant  have  contributed  significantly  to  these  research 
efforts.  In  addition  to  summer  support  for  the  p.i.'s,  the  AFOSR  has  partially  supported 
Basab  Maulik,  a  postdoc  from  Princeton  University;  Karen  Bolinger,  a  1989  Ph.D. 
graduate  under  the  advisement  of  Forest;  Qi  Wang,  a  1991  Ph.D.  graduate  of  Forest's; 
and  Gene  Cao,  a  current  Ph.D.  candidate  of  Forest's  who  will  graduate  in  1992. 

We  summarize  the  main  results  which  have  been  attained  during  this  granting 
period  and  the  publications  which  have  resulted  that  acknowledge  this  AFOSR  grant. 
The  results  will  be  separated  into  an  analytical  category,  a  numerical  category,  and  an 
experimental  category. 

I.  Analytical  Studies 

We  now  have  produced  two  asymptotic  theories  (pointwise  and  averaged)  for 
slender  free  surface  non-Newtonian  jets.  This  analysis:  (a)  reduces  spatial  dimensions 
from  three  to  one;  (b)  yields  a  seven  parameter  family  of  closed  systems  of  pde's  in  one 
space,  one  time  that  govern  slender  free  non-Newtonian  jets;  (c)  establishes  formal 
consistency  to  all  orders  of  perturbation  for  these  1-D  models;  (d)  produces  explicit 
models  of  a  catastrophic  change-of-type  from  hyperbolic  to  mixed  hyperbolic/elliptic 
type.  These  results  are  contained  in  references  [1-6]. 

These  studies  have  spawned  several  new  ongoing  projects  related  to:  numerical 
loss-of-evolution,  change-of-type  in  the  full  3-D  Maxwell  models,  physical  validity  and 
consequences  of  a  mathematical  change-of-type,  numerical  detection  and  consequences 
of  change-of-type,  physically  meaningful  regularizations  that  produce  well-posedness  in 
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an  ill-posed  regime,  the  analysis  of  convergent  and  stable  numerical  codes,  and  molecular 
dynamics  simulations  of  dilute  polymer  solutions.  These  current  projects  in  olve 
collaborations  with  M.  Gunzburger,  VPI,  James  Turner,  Hampton  University,  and  C. 
Jayaprakash,  Department  of  Physics,  Ohio  State  University,  as  well  as  two  Ph.D.  students 
of  Forest,  Q.  Wang  and  J.  Cao. 

n.  Numerical  Studies 

We  have  several  interconnected  numerical  studies,  involving  the  p.i.’s,  Maulik, 
Wang,  Cao,  and  more  recently  Gunzburger  and  Turner. 

Wang  has  focussed  on  the  numerical  simulation  of  specific  asymptotic  model 
systems.  These  codes  compute  well-posed  boundary  value  problems  for  steady  solutions 
based  on  the  characteristics  of  the  pde  systems.  Steady  states  are  computed,  linearized 
stability  of  steady  states  and  effects  of  boundary  conditions  are  computed,  and  the  full 
time-dependent  codes  are  developed  and  run  in  a  neighborhood  of  steady  states  to 
explore  basins  of  stability  and  robustness  with  respect  to  parameters.  The  systems  of 
pde’s  have  change-of-type  behavior,  and  there  are  nontrivial  coding  issues  related  to  this. 
It  is  here  that  mathematically  deduced  physical  regularization  of  a  change-of-type  is 
critical.  We  know  that  in  our  1-D  model  problems  the  retardation  effects  regularize  the 
elliptic  change-of-type  from  high  Weissenberg  numbers,  but  only  in  a  neighborhood  of 
constant  solutions.  The  development  of  a  stable  code  for  the  regularized  equations  is  in 
progress.  This  project  is  prototypical  of  the  3-D  loss-of-evolution  pathologies. 

Cao  has  developed  numerical  codes  which  allow  for  the  investigation  of  the 
asymptotic  validity  of  leading  order  slender  jet  solutions.  Many  solutions  to  "thin 
filament"  models  have  been  presented  in  the  engineering  and  polymer  processing 
literature  since  1969.  Our  work  was  the  first  to  demonstrate  that  these  models  can  be 
obtained  as  the  leading  order  problems  in  a  formally  valid  slender  jet  expansion  (see  §  I); 
now  in  the  paper  with  Cao  [8]  we  are  the  first  to  investigate,  on  a  solution-by-soiution 
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basis,  which  solutions  are  indeed  asymptotically  valid.  Briefly,  we  And  that,  although  all 
of  the  models  are  formally  valid,  particular  solutions  to  these  models  may  be  valid 
(corrections  to  the  leading  order  behavior  remain  small)  or  invalid  (corrections  become 
sufficiently  large  to  disrupt  the  asymptotic  ordering  on  which  the  theory  is  based). 

Maulik  has  numerically  studied  laminar  transition  in  2-D  and  3-D  Maxwell  fluids. 
These  studies  involve  massive  codes  and  the  results  are  now  beginning  to  get  published. 
We  anticipate  comparison  of  these  macroscopic  simulations  with  molecular  dynamics 
simulations  of  dilute  polymer  solutions  by  C.  Jayaprakash,  Department  of  Physics, 
O.S.U.  The  goal  of  this  project  is  to  model  realistic  dynamics  of  the  large  polymer 
molecules  and  ultimately  infer  corrections  to  Maxwell  fluid  models  which  are  connected 
to  realistic  polymer  dynamics. 

III.  Experimental  Studies 

Experiments  have  been  performed  involving  flee  surface  jets  of  four  very 
different  types  of  fluids:  (1)  water,  (2)  viscoelastic  aqueous  biopolymer  solutions  in 
three  different  concentrations  (Xanthan  gum),  (3)  a  viscoelastic  aqueous  surfactant 
solution  (Ethoquad/sodium  salicylate),  and  (4)  water/detergent  solutions. 

The  experiments  imply  that: 

1)  robust,  slender  jet  profiles  exist,  so  that  the  assumptions  of  our  asymptotics  are 
supported; 

2)  noticeable  effects  due  to  material  (viscoelastic)  properties  are  documented  in 
otherwise  identical  experiments. 

An  open  problem  in  rheology:  How  does  one  measure  relaxation  time  Aj  in 
eiongational  flows?  (Measurement  of  elongational  response  in  fluids  has  proven  difficult 
because  regions  of  elongational  flow  are  in  general  difficult  to  isolate  and  involve  free 
surfaces,  which  complicates  the  analysis.  Elongational  response  cannot  be  inferred  from 
the  more  easily  measured  shear  response  of  a  fluid.)  Our  asymptotic  models  for  free  jet 
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flow  (which  is  elongational)  predict  experimentally  measurable  effects  as  functions  of  A  j 
for  other  parameters  fixed.  This  poses  an  intriguing  inverse  problem  from  our 
analysis:  we  infer  Aj  by  comparing  experimental  properties  of  the  free  jet  flow  (e.g.,  the 
wavelength  and  amplitude  of  the  observed  nonlinear  spatial  oscillation  in  the  steady 
profile  of  a  jet  issuing  from  an  elliptical  orifice,  and  the  growth  or  decay  of  these 
quantities  in  successive  cycles  down  the  jet)  with  the  corresponding  value  of  Aj  that 
produces  the  same  properties  in  the  1-D  model. 
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